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q . ABSTRACT 

The dwarf nova GW Librae is the first cataclysmic variable discovered to have a primary in a 
r/3 , white dwarf instability strip, making it the first multi-mode, nonradially-pulsating star known 

d • to be accreting. The primaries of CVs, embedded in hot, bright accretion discs, are difficult 

to study directly. Applying the techniques of asteroseismology to GW Librae could therefore 
give us an unprecedented look at a white dwarf that has undergone ~ 10 9 years of accretion. 
1 However, an accreting white dwarf may have characteristics sufficiently different from those 

J-j | of single pulsating white dwarfs to render the standard models of white dwarf pulsations in- 

valid for its study. This paper presents amplitude spectra of GW Lib from a series of observing 
campaigns conducted during 1997, 1998 and 2001. We find that the dominant pulsation modes 
cluster at periods near 650, 370 and 230 s, which also appear in linear combinations with each 
other. GW Lib's pulsation spectrum is highly unstable on time-scales of months, and exhibits 
clusters of signals very closely spaced in frequency, with separations on the order of 1 filiz. 
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1 INTRODUCTION 

The dwarf nova (DN) GW Librae was the first cataclysmic vari- 
able (CV) discovered to contain a nonradially pulsating white dwarf 
(DAV, or ZZ Ceti star; van Zyl et al. 2000). The existence of a DAV 
within an interacting binary system has some potentially signifi- 
cant implications for the study of accreting stars. White dwarfs in 
CV systems, with their long accretion histories and many novae 
eruptions during their ~ 10 9 -10 10 -year lifetimes, may have very 
different internal structures from single white dwarfs (Townsley 
& Bildsten 2002). However, the primaries of CVs are difficult to 
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study: the flux from CVs tends to be dominated by accretion pro- 
cesses. Nonradial pulsations, which carry a wealth of information 
about stellar structure, composition and evolution (e.g., Winget et 
al. 1991, 1994), therefore potentially provide us with a unique and 
very powerful tool with which to study CV evolution and the ac- 
cretion process. 

In addition, the existence of an accreting DAV also has some 
exciting implications for the study of single white dwarfs. As DA 
white dwarfs cool on evolutionary time-scales, taking 5-10 x 10 s 
years to pass through the DAV instability strip (Wood 1992), their 
pulsations change in ways not yet fully understood, becoming more 
unstable and prone to nonlinear effects. In the case of a DAV in a 
CV system which undergoes periodic DN outbursts, it may be pos- 



2 L. van Zyl et al. 





* » 


* f #v 






0.23 0.^4 






0.28 


0.3 












0.32 0.34 


— i- 


— 1 ' ' r- 

0.36 


0.3E3 C 


4 


*A* ^ ¥ f V% 






x **** x 





0.42 0.44 0.46 48 0.5 

IIJD - 2450965 



Figure 1. A sample of GW Lib's light-curve, obtained with the SAAO 0.75- 
m telescope on 31/05/1998. 

sible to observe the star move through the traditional DAV insta- 
bility strip in a matter of only months or years, ~ 10 8 times more 
rapidly than single DAVs, as the white dwarf photosphere cools 
between DN outbursts. 

1.1 The first accreting multi-mode pulsating star 

GW Librae is a faint DN which called attention to itself for the first 
time in 1983 when its brightness increased by ~ 9 mag (Gonzalez 
1983) - its only observed outburst. The large amplitude of the out- 
burst initially led to GW Lib being misclassified as a nova; how- 
ever, subsequent spectroscopy (Duerbeck & Sitter 1987; Ringwald, 
Nay lor & Mukai 1996) showed it to be a system with a very low 
mass-transfer rate. GW Lib is almost certainly a member of the 
class of long outburst-interval, low mass-transfer-rate DNe known 
as the WZ Sagittae stars. 

GW Lib's faintness (V~18) precluded interest in it for 14 
years following its outburst, until we chanced upon its remark- 
able properties while conducting a high-speed-photometry survey 
of faint Southern Hemisphere CVs. Its spectrum and its 1983 WZ- 
Sge-style superoutburst show beyond doubt that GW Lib is a CV, 
but the Fourier transform of its light-curve (Fig. 1) resembled that 
of a nonradially pulsating single white dwarf. Nonradial pulsations 
had never been observed in an accreting white dwarf before (ex- 
cept possibly a single mode in the DOV primary of the CV AM 
CVn [Solheim et al. 1998]); it had been assumed either that ac- 
cretion would keep the primaries of CVs too hot to pulsate, or 
that He in the accreted material would suppress pulsations. The 
DAV (or ZZ Cet) instability strip for (non-accreting) hydrogen- 
atmosphere white dwarfs (DAs) occurs between T e f?s of ~ 11,000 
and ~ 12,500 K (Koester & Holberg 2001), depending on the mix- 
ing length prescription used to describe convection. The coolest DN 
primaries have surface temperatures of typically 15 000 K (Sion, 
Urban & Lyons 2001). 

A spectroscopic study of GW Lib by Thorstensen et al. (2002) 
found an orbital period of 76.78 min - the shortest known orbital 
period for a CV with a hydrogen-rich donor (V485 Cen has a pe- 
riod of 59 min, but its companion has a low hydrogen content [Au- 
gusteijn et al. 1996]). This orbital period supports GW Lib's status 
as a WZ Sge star. The WZ Sge stars, or 'TOADs', are the old- 
est DNe, and are characterized by very low mass-transfer rates, 
long outburst recurrent times and short orbital periods (e.g. Bai- 
ley 1979; O'Donoghue et al. 1991; Howell, Szkody & Cannizzo 
1995). Along with the magnetic AM Her CVs, which also have 



low mass-transfer rates, WZ Sge stars are therefore the most likely 
CVs to harbour DAVs, because of their age (white dwarfs take 0.5- 
5 Gyr to cool to 12,500 K, depending on their mass [Chabrier et al. 
2000]), and because of their low mass-transfer rate (so that com- 
pressional heating of the white dwarf's cores through accretion is 
small [Townsley & Bildsten 2001]). 

Despite extensive photometry, no orbital modulation has been 
detected in GW Lib's light-curves, suggesting that the system has 
a low inclination. From the width of GW Lib's emission lines, 
Thorstensen et al. (2002) infer an inclination of ~ 1 1°. Thorstensen 
et al. also find an absolute magnitude for GW Lib of Mv ~ 11.5, 
a distance of ~ 125 pc and a significant proper motion of 66 ± 
12 mas y _1 . 

Szkody, Desai & Hoard (2000) found that spectra of GW 
Lib were roughly consistent with an effective temperature of 
11,000± 1,000 K (assuming that 100% of the flux came from the 
white dwarf), which put GW Lib in the DAV instability strip. 
However, Szkody et al. (2002) then obtained a more accurate 
white dwarf surface temperature, with HST UV spectroscopy, of 
14,700 K, which puts it well outside the traditional DAV instabil- 
ity strip. This may imply a possible temporary heating of the white 
dwarf atmosphere due to an unobserved outburst (while the core 
and mean temperatures are much lower), or that the instability strip 
for accreting DAVs is different from that of single DAVs. Indeed, 
the notion of an 'instability strip' for accreting white dwarfs may 
be meaningless: instability strips are usually associated with groups 
of stars moving by single star evolution through a T e fi interval, but 
in accreting systems the mass transfer onto the white dwarf would 
very likely introduce very important additional parameters. 

GW Lib's pulsation spectrum shows the unstable behaviour 
typical of the cool, large amplitude DA white dwarf pulsators. As 
white dwarfs evolve from the hot to the cool end of the DAV in- 
stability strip, their pulsation spectra become increasingly more 
complex and unstable, the amplitudes of their pulsation modes of- 
ten changing dramatically from month to month (Kleinman et al. 
1998). The amplitudes of the oscillations in GW Lib's light-curve 
are very small (~ 5-10 mmag). However, as we expect approxi- 
mately half the light from the system to be coming from the accre- 
tion disc (based on a comparison with the DN Z Cha), the intrinsic 
amplitude of the pulsations is likely to be greater. 

1.2 Problems specific to accreting pulsators 

As a result of its birth in a common envelope and its subsequent 
long history of accretion, it would not be surprising if GW Lib's pri- 
mary were to have characteristics, and therefore a pulsation spec- 
trum, different from those of single white dwarfs. We discuss below 
some issues which we believe may be of significance in the case of 
an accreting DAV. 

1.2. 1 Accretion-induced spherical asymmetries 

If a CV primary has a weak or absent magnetic field (as is the 
case of GW Lib, which does not appear as an X-ray source in 
the ROSAT All-Sky Survey Faint Source Catalogue [Voges et al. 
2000]), material from the disc will accrete onto the equatorial re- 
gions of the star. The white dwarf therefore develops an equato- 
rial band which has a different temperature, rotational and chem- 
ical composition structure to the cooler, higher-latitude regions of 
the star. The effect is especially pronounced after superoutbursts. 
Sparks et al. (1993) find that after WZ Sge's superoutbursts, the 



The Pulsating Primary of the CV GW Librae 3 



Period (s) 

2000 1000 700 500 400 300 200 




Frequency (mllz) 



Figure 2. The amplitude spectra of GW Lib from 6 observing campaigns 
from March 1997 to May 2001. The light-curves were prepared by fitting 
and removing linear and low-frequency trends. 

hot accreted material is confined to a broad equatorial band on the 
primary for ~10 years. Nonradial modes are sensitive to such lat- 
itudinal differences within the stellar envelope (e.g. Schou et al. 
1998). However, single DAVs are spherically symmetric to a very 
high degree, and therefore their pulsation spectra have been free of 
these complexities. 

As yet, no-one has tried to model the pulsational behaviour 
of an accreting DAV, so at present we can only speculate as to the 
effects of accretion on its pulsation spectrum. If accretion-induced 
spherical asymmetries extend to the depth of the pulsation driving 
zone, the pulsation eigenfunctions could be dramatically affected. 

In addition, it is possible that an accreting DAV may exhibit 
pulsation modes with / > 2: in single DAVs, geometric cancella- 
tion effects result in modes with I > 2 having amplitudes too low 
to be detected. However, an accreting DAV would have a differ- 
ent surface temperature distribution, which could alter the surface- 
temperature geometry of modes with I > 2 in such a way that they 
may become observable. 

1.2.2 Rapid Rotation 

Single DAVs rotate very slowly, with periods on the order of sev- 
eral hours or longer (e.g. Koester et al. 1998, Bradley 2001), and 
their pulsations have periods on the order of hundreds of seconds. 
Therefore 

where cu {k , m) are the pulsation eigenfrequencies and O, is the stel- 
lar rotation frequency, and the leading effects of rotation on the 



oscillations are linear in fi. The effect of slow rotation is to lift 
the degeneracy of the , m) eigenfrequencies, splitting them into 
multiplets, where m = 0, ±1, ±Z. 

CV primaries, however, typically have rotation periods much 
shorter than those of single white dwarfs, because they are spun 
up by the accretion process (e.g. Warner 1995). The rotation peri- 
ods for primaries in non-magnetic CVs are on the order of 100 s 
(Warner [1995] and references therein). In this case, 



and rotation can no longer be dealt with as a small perturbation, 
and therefore the pulsation modes can no longer be described by a 
single spherical harmonic Y^O, (f>). Coriolis forces become of the 
same order as buoyancy forces which provide the restoring forces 
for p-modes, and therefore become a new restoring force for waves 
of a very different physical nature, which can also be trapped in 
the star and form global modes (Vorontsov 1993). The geometry 
of the Coriolis forces is significantly more complex, making the 
classification of the new modes difficult. 

1.2.3 Non-stationary pulsation spectra 

During a superoutburst, the temperature of a WZ Sge primary's 
photosphere can increase by a factor of 2 or 3, and take ~ 10 years 
to return to its quiescent value (Sion & Szkody 1990; Sparks et 
al. 1992; Long et al. 1994). The characteristics of DAV pulsations 
are sensitive to temperature, and therefore superoutbursts are likely 
to have an effect on the pulsation spectrum of a CV-primary DAV. 
Frequency and amplitude modulation as a result of temperature- 
induced changes in the envelope could make the pulsation spectrum 
very difficult to unravel. 

In addition to steady temperature changes of the photosphere 
resulting from superoutbursts, there may also be shorter-timescale 
accretion-induced temperature perturbations which may affect the 
pulsation modes. Flickering in DNe light-curves indicates that ac- 
cretion is very 'blobby' (Warner 1995). Blobs of material accreting 
onto a star unstable to pulsations may behave like a stochastic ham- 
mer applied to an oscillator, inducing constant phase changes in the 
oscillations, which would also complicate the pulsation spectrum. 



2 THE PULSATION SPECTRUM OF GW LIB 
2.1 The observations 

We conducted six observing campaigns on GW Lib: in March, 
April and September 1997, May and June 1998, and May 2001. 
A summary of the observations is given in Table 1. The bulk 
of the data were obtained from the South African Astronomical 
Observatory (SAAO) in Sutherland, South Africa, using a high- 
speed, frame-transfer CCD photometer purpose-built for observ- 
ing blue, rapidly-varying stars like pulsating white dwarfs and 
CVs (O'Donoghue 1995). A large part of the data we discuss has 
appeared in the proceedings of the Fifth Whole Earth Telescope 
Workshop (van Zyl et al. 2000). 

In order to attempt to resolve GW Lib's pulsation spectrum, 
we conducted a multi-longitude observing campaign in May 1998. 
The participating observatories were Mt John Observatory (New 
Zealand), Mt Stromlo & Siding Spring Observatory (Australia), 
CTIO (Chile) and SAAO. All sites made use of CCDs, and ob- 
servations were made in white light (due to GW Lib's faintness), 
but only SAAO was able to use a blue-sensitive, back-illuminated 
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Figure 3. Expanded plots of the 650-s, 370-s and 230-s regions of GW Lib's amplitude spectra. The fourth column gives the spectral windows for each 
observing campaign. The unstable nature of GW Lib's pulsations is clearly apparent. Our best data, with the highest signal-to-noise, the longest run-lengths 
and the cleanest spectral window, are the May 1998 observations. However, in May 1998 GW Lib's pulsation spectrum, apart from a new mode near 230 s, 
showed fewer active modes than in April 1997. 



chip. Very poor weather at all the sites apart from SAAO resulted 
in us obtaining only two half-nights of data from Mt John, only one 
night from Mt Stromlo, and no usable observations from CTIO. 
From SAAO we were able to observe on 13 of 14 nights, and on 
most of these nights were able to obtain good quality, long light- 
curves. 

In order to extend the baseline of the May 1998 dataset (and 
thereby increase the resolution of its Fourier transform), we ob- 
tained 3 good runs at SAAO on 3 successive nights in June 1998, 
separated from the May 1998 data by 14 days. We also obtained 
short light-curves from an additional 3 nights in June 1998 from 
McDonald Observatory (USA), 9 days after the end of the June 



SAAO run. However, adding these data to the SAAO dataset com- 
plicated the spectral window to an extent which negated any ben- 
efit from further extending the baseline. Therefore, in what fol- 
lows, 'May98' refers to the May 1998 multisite campaign, and 
'May&June98' refers to the May98 dataset with its baseline ex- 
tended with the inclusion of the 3 nights of SAAO data obtained 14 
days later. 

2.2 The amplitude spectra 

The amplitude spectra presented in this paper are Fourier trans- 
forms (FTs) of GW Lib's light-curves, obtained using an algorithm 
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SAAO = South African Astronomical Observatory, MTJ = Mt 
John Observatory (New Zealand), MSSSO = Mt Stromlo & Sid- 
ing Spring Observatory (Australia), MCD = McDonald Obser- 
vatory (USA). 

The May 1998 observations are our best data, with the longest 
run-lengths, highest signal-to-noise and cleanest spectral win- 
dows. 

These observations from June 1998 are our next best dataset, 
and when combined with the May 1998 dataset to extend its 
baseline, make up the 'May&June98' dataset referred to in the 
text. 
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Figure 4. The 230-s mode. In the amplitude spectrum of the May98 data 
alone, this feature appears to be a single mode (left). However, with the 
addition of 3 nights in June 1998, the baseline is sufficiently extended to 
reveal a second, low-amplitude mode after prewhitening with the dominant 
mode (right). (The 'bump' left of centre is a linear combination mode.) 



for calculating discrete FTs of unequally-spaced data proposed by 
Deeming (1975), and with a short-cut proposed by Kurtz (1985). 
In Fig. 2 we present amplitude spectra of GW Lib from each of the 
six observing campaigns. Fig. 3 shows the regions of interest in the 
amplitude spectra in greater detail. 

GW Lib's pulsation spectrum is at its most active in April 
1997, when there appear to be two modes active near 650 s and 
three near 370 s (Figs. 2 and fig:details). The structure of the modes 
in the April 1997 Fourier transform is a great deal more complex 
than the window function (Fig. 3), indicating that each mode con- 
sists of more than one component. However, the April 1997 dataset 
does not have the baseline, nor the duty cycle, necessary to resolve 
the fine structure. 

By the time of the May98 multisite campaign, however, GW 
Lib's pulsation spectrum had changed dramatically. Three of the 
five pulsation modes present near 650 s and 370 s in April 1997 
had disappeared by May 1998, and a new one near 230 s, not previ- 
ously seen, had become detectable. For convenience, we shall refer 
to the three principal structures in the May98 and May&June98 
amplitude spectrum, at 650 s, 370 s and 230 s, as u) 650 , ui 370 and 
cj 230 . However, we show below that each of these is a multiplet 
consisting of more than one frequency. 

In Table 2 we list all the signals present in each amplitude 
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Figure 5. O-C diagrams (left) and amplitudes (right) of the three princi- 
pal structures,w 230 , u> 370 and lo 6so , in the May&June98 amplitude spec- 
trum. The points give the O-C and amplitude of the structures for the first 
and second half of each night, calculated using the mean frequency of each 
structure's components from Table 2. 



spectrum. Each signal consists of a set of aliases (ambiguities 
caused by gaps in the data), and Table 2 lists only the highest- 
amplitude alias for each signal. There is ambiguity in the aliases 
(the highest-amplitude alias is not necessarily the true frequency), 
so each frequency uj in Table 2 should be read as a set of frequen- 
cies co + n x 11.57/iHz, where n = ±1, ±2, ±3, and 1 1.57 /iHz 
is the dominant 1-day alias. 

The procedure for determining the amplitudes and frequen- 
cies of the signals in the amplitude spectrum is as follows. First, 
the mean and second-order trends (due to sky transparency) are re- 
moved from the light-curves, and then the amplitude spectrum is 
calculated. Then, assuming there are no aliasing ambiguities, the 
frequency of the highest peak is fitted to the data and subtracted (a 
process known as 'pre-whitening'). The amplitude spectrum of the 
residuals is then computed and the process repeated until the data 
are consistent with noise. The frequencies found in this way (lin- 
ear combinations as well as principal frequencies) are then fitted 
simultaneously to the full data set by a non-linear least-squares fit, 
which produces our final set of frequencies and amplitudes. 



2.3 The mode at 230 s 

Ironically, the simplest feature of the May98 amplitude spectrum is 
also the most obviously unstable: there is no trace of the mode near 
230 s in any of the 1997 amplitude spectra, although it is strong in 
2001. At first glance, u) 230 (Fig. 4, upper panels) appears identical 
to the spectral window in both the May98 and the May&June98 
FTs. The lower panels of Fig. 4 show the FTs after prewhitening 
with this signal. In the May&June98 spectrum, with its extended 
baseline, an additional, low-amplitude component is revealed. 

If the 230-s signal is a rotationally-split doublet (the most 
common reason for splitting of DAV modes), the frequency dif- 
ference between the components, 0.79 /iHz, is remarkably small. If 
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Figure 7. Linear combination modes in the residuals of the May&June98 
light-curves after prewhitening with all the components of the 650-s, 370-s 
and 230-s modes. 



this is an I — 1 mode (the most common modes in DAVs), the im- 
plied rotation period of the star is 7.3 days (from the properties of 
pulsational g-modes; see for example Winget et al. [1991,1994]). 
This rotation period is much longer than what we would expect for 
a CV primary. 

A splitting of 0.79 /iHz corresponds to a beat period of ~ 15 
days, which is very close to both the baseline of the May98 dataset, 
and the gap between the May and June data. However, there are 
no splittings equal to 0.79 /iHz in the spectral window. Further ev- 
idence that w 230 is a closely-spaced doublet comes from changes 
in its O-C phases and amplitude over the observing campaign. 
In Fig. 5 we show the O-C and amplitudes of the three principal 
regions of power, o> 230 , uj 370 and w 6S0 , over the baseline of the 
May&June98 dataset. The points were calculated using the domi- 
nant signal in u> 230 and an average of the frequencies for each of 
u> 370 and lj 650 listed in Table 2. From the top panels in Fig. 5 it 
can be seen that the O-C curve of uj 230 is not flat (Fig. 5), which 
could indicate that w 230 consists of two components. Fig. 5 also 
shows that the amplitude of a; 230 appears to be varying sinusoidally, 
which could indicate the existence of two closely-spaced compo- 
nents beating against each other. 



2.4 The multiplets near 370 s and 650 s 

The other principal regions of power in the pulsation spectrum of 
the May98 dataset - the structures near 650 s and 370 s - are a 
great deal more complex than co 23a . Fig. 6 shows the 370- and 650-s 
regions of the May98 and May&June98 amplitude spectra in detail. 

ui 370 appears to consist of two triplets, each with compo- 
nents spaced ~ 1 pHz apart (Table 2). One triplet, centred at 
2645.90 /iHz or 377.94 s (indicated by the open arrows in Fig. 6, 
which show the three highest-amplitude aliases of the complete 
triplet), has three closely-spaced components well resolved in both 
the May98 and May&June98 FTs. No feature of the spectral win- 
dow corresponds to the frequency spacing of the triplet compo- 
nents, so we do not believe that they are an observational artifact. 
u> 370 's second structure appears to be an unresolved triplet centred 
at 2662.94 /iHz or 375.52 s (indicated by the solid arrows in Fig. 6). 
Prewhitening this structure produces 3 signals also ~ 1 /iHz apart. 
We fitted these three signals simultaneously to the May98 dataset, 
the May&June98 dataset, and two additional datasets containing 
only the first or the second half of each night. In all four cases, 
this feature was consistent with a triplet with components ~ 1 /iHz 
apart. 

The presence of both a resolved triplet and an unresolved 
triplet whose components have the same spacing is not inconsis- 
tent: if a light-curve is full of substantial gaps (due to daylight or 
bad weather, etc), then the resolvability of closely-spaced signals 



Table 2. The major periodicities (excluding linear combinations) detected in GW Lib's amplitude spectra from 1997 to 2001. Low-frequency signals have been included, but the list is 
likely to be incomplete as low-frequency trends were removed from the light-curves in order to bring down the level of the noise. 
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Figure 6. Details of the amplitude spectra of the May98 (upper panels) and May&June98 datasets (lower panels). The power in these regions is dominated by 
very line frequency splitting, on the order of ~ 1 £tHz, not all of which can be explained in terms of artifacts generated by the window function. The 'mode' 
we call u> 37a (left-hand panels) appears to consist of two triplets, each with components spaced ~ 1 £tHz apart. One triplet (the three highest-amplitude aliases 
indicated by the open arrows) is centred at 2645.90 /iHz or 377.94 s. The second structure (the three highest-amplitude aliases indicated by the solid arrows) 
appears to be an unresolved triplet centred at 2662.94 ^Hz or 375.52 s. 



whose beat period is approximately equal to the total run length 
(as is the case here) depends on how the data and the gaps are dis- 
tributed across their beat cycle. If the data cover regions where the 
signals are out of phase with each other, the FT is far more likely 
to resolve the individual signals than if the signals are in phase. We 
have modelled pairs of triplets made up of closely-spaced sinusoids 
sampled in the same way as our data, and we are able to reproduce 
the appearance of w 370 very well. 

w 650 , on the other hand, is a partially-resolved tangle of at 
least 7 closely- but unevenly-spaced signals (Table 2). It is diffi- 
cult to know whether these signals represent individual pulsation 
modes active within the white dwarf, or whether we are seeing 
the effects of non-stationary behaviour in the pulsation spectrum 
(i.e., frequency, amplitude or phase modulation). The O-C curves in 
Fig. 5 (middle and bottom panels) show noisy modulations, which 
support our findings that lo 370 and u) g50 appear to be made up of 
several different signals. 



2.5 Linear combinations 

Linear combination modes are frequently observed in the pulsation 
spectra of DAVs. They arise as a result of non-linear mixing of the 
real eigenmodes: the thickness of the convection zone is modulated 



by the pulsations and the signal is not only diminished and delayed, 
but also distorted. This latter distortion gives rise to the power seen 
at linear combinations (Brickhill 1992). 

We have detected four clusters of signals at frequencies cor- 
responding to linear combinations in the May98 and May&June98 
amplitude spectra. In Fig. 7 we plot the FT of the residuals of the 
May&June98 dataset after prewhitening with lj 650 , oj 370 and ui 230 . 
The amplitudes of these linear combination modes are very low, 
only a factor of two above the noise level. If they are present in the 
1997 and 2001 amplitude spectra, which have greater noise levels, 
they would therefore not be detectable. 

Three of the four linear combination regions lie very close 
to the dominant signals: the combination (uJ 23a — w 370 ) is only 
^40 /iHz from uj 6S0 , (uj 230 — io eso ) is ~30 /iHz from co 370 , and 
(^650 + w 37o) is ~40 nHz from lo 250 . The proximity of the linear 
combinations to the dominant signals is unfortunate; at these small 
separations the window functions partially overlap, which almost 
certainly contributes to the confusion in lo 650 and co 37a . 

Our choice of whether a signal is due to a real or combina- 
tion mode is dictated by their relative amplitudes: the combination 
mode is defined to have a smaller amplitude than the constituent 
real modes. Despite their low signal-to-noise, we discovered 24 sig- 
nals in the linear combination regions that corresponded to better 
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Table 3. Linear combinations. 



Combination Calculated Observed Auj 
GuHz) QiHz) QiHz) 

0.09 
0.46 
0.07 
0.81 
0.41 
0.26 
2.02 
2.60 
1.20 

0.04 
0.06 
0.41 
0.66 
0.41 
0.53 
0.48 
0.13 
0.62 
1.21 

0.26 
0.43 
0.26 
1.61 
0.14 

1.00 
0.08 
0.10 
1.10 
0.56 
0.40 



cii = the 6 highest-amplitude components at 650 s 
bi = the 6 highest-amplitude components at 370 s 
c = the dominant component of the 230s doublet 
Aw = observed - calculated 



than 1 /jHz (some to better than 0. 1 /iHz) to the set of 43 combina- 
tions calculated from the highest-amplitude components of uj B50 , 
u> 370 and ui 23a (Table 3). By calculating the fraction of the fre- 
quency axis that is filled by possible linear combinations, we esti- 
mate that only 6 coincidences would be expected by chance align- 
ments. Thus, we have confidence that the fine structure in lo 650 and 
w 370 is real: a signal is more likely to be real if it appears in linear 
combinations as well (Kleinman et al. 1998). 



3 DISCUSSION 

GW Lib appears similar to the cool DAVs (its pulsation periods 
are consistent with those expected for non-radial g-modes in single 
DAVs), which have highly unstable pulsation spectra, with modes 
disappearing and reappearing on monthly or yearly time-scales, 
and which therefore require several years - sometimes decades - 
of multisite observing campaigns in order that enough modes are 
detected for mode identification (and therefore asteroseismology) 
to be possible. Many more modes will need to be identified before 



we can begin to unravel GW Lib's pulsation spectrum. But is there 
anything we can say about its pulsation modes at this early stage? 

Already, there are hints that GW Lib's pulsation spectrum is 
different from those of single DAVs. The frequency splitting is 
much finer than the multiplet splitting observed in single DAVs, 
and if caused by rotational splitting, would imply that GW Lib's 
primary is rotating an order of magnitude more slowly than sin- 
gle DAVs, which for a CV primary is improbable. In addition, if 
we consider cu 37Q 's two triplets to be modes of same / (because the 
spacing of their components is the same) then the separation be- 
tween modes is too small. From Bradley's (1996) models, the min- 
imum average separation for DAV modes of I = 1 is ~ 34 s, and 
~ 20 s for I = 2 modes, while the two triplets within io 37a are sep- 
arated by only ~2 s. Two triplets of different I could have a small 
separation, but then their component spacings would be different. 

The April 1997 amplitude spectrum (Fig. 3) shows (at least) 
three 'modes' near 370 s (each mode appears to consist of more 
than one structure, so may therefore be composed of two or more 
multiplets, like ui 370 ). The right-hand April97 mode (at 378.9) dif- 
fers by only ~ 6 /iHz from the left-hand triplet in lo 370 , so - given 
the poorer resolution of the April97 dataset - these may be the same 
mode. This still gives us four, perhaps more, modes lying within a 
20-s interval. Even if we imagine them to be modes of alternating 
I = 1 and / = 2, they are still too close together. 

These awkward period spacings, both the fine splitting within 
apparent modes and the mode-to-mode spacing, show that it is go- 
ing to be difficult to understand GW Lib's pulsation spectrum using 
the framework appropriate to single, non-accreting DAVs. It may be 
that the assumptions implicit in the above paragraphs (i.e., that the 
modes are all of high radial order, and that only modes of I = 1 and 
I = 2 are present) don't hold for accreting DAVs. It may be that in 
GW Lib we are observing, for the first time in a DAV, modes of 
I > 2. In §1.2.2 we point out that geometric cancellation may work 
differently for stars with hot equatorial bands, allowing modes of 
higher I to be visible. A possible way to identify whether GW Lib's 
modes are of I — 3 or higher would be to obtain simultaneous 
UV and optical light-curves; higher-order modes have higher UV 
amplitudes because limb-darkening helps to negate the effects of 
geometric cancellation at these wavelengths (Robinson et al. 1995; 
Kepler et al. 2000). In the case of GW Lib an appropriate model 
for the fractional pulsation amplitudes would have to be developed, 
taking into account possible accretion-related differences in sur- 
face temperature distribution and their effect on limb-darkening 
and geometric cancellation, before this mode-identification tech- 
nique could be applied. 

If we are to understand GW Lib's pulsation spectrum, many 
more pulsation modes will need to be identified, which will require 
a substantial observational effort. In addition, the behaviour of g- 
modes in accreting, rapidly rotating white dwarfs with latitude- 
dependent temperature, angular momentum and chemical compo- 
sition structures in their outer layers will need to be modelled, as 
will the effects of DN outbursts and superoutbursts on the stars' 
eigenfrequencies. 
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